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Abstract 

Feasibility of the measurement of the trilinear self-couplings of the Higgs boson is 
studied. Such a measurement would experimentally determine the structure of the Higgs 
potential. Full hadronic and semi-leptonic final states of the double-Higgs strahlung have 
been investigated. 



1 Introduction 



In the framework of the standard model, the generation of mass occurs through the Higgs mech- 
anism. This mechanism assumes a Higgs potential, V($), which behaves as V(<&)=\(<$> 2 -^v 2 ) 2 , 
where <fi is an isodoublet scalar field, and t>~246 GeV is the vacuum expectation value of its 
neutral component. A determination of the Higgs boson mass, which satisfies m 2 1 =4Av 2 at tree 
level in the standard model, will provide an indirect information about the Higgs potential 
and its self-coupling A. The measurement of the trilinear self-coupling Aj 1 j 1 j 1 =^ \v offers an 
independent determination of the Higgs potential shape. In case of disagreement, the compar- 
ison between these two measurements would give an interesting clue about new physics, like 
super-symmetry for instance. Anyway, the reconstruction of the Higgs potential is an essential 
step in the experimental validation of the Higgs mechanism that lies at the core of the standard 
model. 

The aim of our study is to demonstrate how to measure the trilinear self-coupling A^j^ in 
e + e~ collisions at center-of-mass energies delivered by the future Linear Collider. The note is 
organized as follows : the framework of the study is defined in Section |2|. The Monte Carlo 
simulation is reported in Section |3|. Analyses and results are described in Section f| and [|. 

2 Signal 

The trilinear Higgs self-coupling, Av^h, could be extracted from the measurement of the cross- 
section of each of the following processes : double Higgs-strahlung (e + e~ — > Zhh) or WW 
double-Higgs fusion (e + e~ — > z/ e z/ e hh) [|TJ. Figure [I] indicates the major Feynman diagrams 
involved in the first process. 




Figure 1: Feynman diagrams involved in the e + e — >hhZ cross-section via the double Higgs-strahlung. 

For light Higgs boson masses, the Higgs boson decays predominantly in a bb pair. When 
the Z decays in a lepton pair, the final state is hhZ^bbbM + £~. With a high b content and two 
leptons, this topology produces an easy signature but represents only ~ 8% of the total final 
state. By contrast, the hhZ— »bbbbqq final state benefits from a high statistics with ~ 60% of 
the final states but requires a more complicated analysis. 

In the minimal super-symmetric extension (MSSM) of the standard model, the Higgs sector 
is richer with three neutral Higgs bosons (h, H and A). The Higgs potential is therefore more 
intricate and further trilinear self couplings are possible beyond A^j^ : A^g, Agj^jj, Aj^A' 
and Aj^A- At the tree level, these depend on two parameters (tan/3 and for instance). 
The radiative corrections are driven mainly by the top quark, so that the stop mass effectively 
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constitutes a third one. Studies performed in the parameter space of the MSSM are described 
in Ref . [l|, [J . Here, let us simply recall on Fig. |] the sensitivity of A^j^ as a function of with 
different m~ t assumptions. For fixed m^, the value of A^j^ reaches its standard model value 
for large (decoupling limit), and the dependence on only shows up for smaller values of 
7T7.A- However, the topologies obtained will be similar to the standard model ones, except for 
particular parameter choices where the Higgs boson decay into a bb pair is suppressed. Thus, 
the experimental study we develop in this note, though a priori restricted to the standard model 
framework, is also relevant for the MSSM. 




Figure 2: MSSM trilinear Higgs self-coupling (A^h) as a function of the mass of the pseudo- 
scalar Higgs boson mass (mj^) and different assumptions of mi from 300 GeV/c 2 to 1 TeV/c 2 with 
a 100 GeV/c 2 step size; the horizontal line corresponds to the standard model value of A^hh- -A 
100 GeV/c 2 Higgs boson mass has been assumed. 



If the mass of the Higgs boson is assumed to be light (~ 120 GeV/c 2 ) then the best center- 
of-mass energy is around 500 GeV [0 and the double Higgs-strahlung process cross-section is 
one order of magnitude greater than the WW fusion one. The study presented here will be 
performed for m h =120 GeV/c 2 and -^=500 GeV and it is restricted to the standard model 
framework. 

Nevertheless, with mh=120 GeV/c 2 and a/s=500 GeV the e + e~— >hhZ cross-section is still 
very tiny with a 0.18 fb value. Namely, for an expected integrated luminosity of 500 fb _1 , only 
93 signal events are produced. The difficulty is increased by the the huge background sources, 
as detailed in next Section. 
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Table 1: Cross-sections for signal and background processes, Monte Carlo statistics and simulated 
luminosity (C s im)- 

3 Simulation 

The signal event samples have been simulated with the GRACE generator [Q. The background 
sources are either two fermions final state processes as e + e~— >Z7, or four fermions processes 
like e + e~— >W + W~, ZZ, We^, Zee, hZ. The latter has been simulated with HZHA generator [|| 
while PYTHIA || was used for the others processes. Six fermions final state processes like 
e + e~— >WWZ or ZZZ are also background sources and GRACE was used to simulate them. Ta- 
ble [l] summaries the size of the samples used and the cross-section of each process. Since 
e + e~— ►tt and e + e~— >Wtb will yield major contributions to the background, an emphasis has 
been performed and a simulated luminosity greater than 2000fb _1 was generated for such pro- 
cesses as well as for six fermions processes while for the others a simulated luminosity greater 
than 500fb _1 was used as reported in Table [l|. Even if the cross-section is small, e + e~— >tth 
process is also a potential source of background and this process was considered. 

The detector simulation was performed with a Parametric Monte Carlo J7J. The tracking 
system, immersed in a 3 tesla magnetic field, consists on a vertex detector (VDET) surrounding 
the beam tube at ~ one centimeter radius followed by a time projection chamber (TPC). The 
VDET impact parameter resolution is assumed to be ~5/im in both rz and r0 views. Only 
charged particles with a transverse momentum greater than 0.2 GeV/c are reconstructed. The 
tracking efficiency is assumed to be 99.7%. The tracking system is complemented with a forward 
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tracker and muon chambers. The b-tagging is performed through a parametrisation derived 
from full reconstruction ||. 

At low angle, a luminometer (LCAL) enforces the hermeticity of the detector. Energy 
and direction of the photons and neutral hadrons are measured thanks to the electromagnetic 
(ECAL) and hadronic (HCAL) calorimeters. The intrinsic energy resolutions taken into account 

are AE / £=10.2%/ \j ' E{GeV) and AE/E=A0.5%/^/E(GeV) respectively for ECAL/LCAL and 
HCAL. The energy threshold below which no neutral particles is reconstructed in calorimeters 
are 200 MeV and 500 MeV for ECAL and HCAL respectively. Following studies performed 
with full simulation of high granularity calorimeters ||, a jet energy resolution of 40%/ \J Ej et 
could be achieved. Such a resolution has been assumed in the analysis described hereafter. 

4 Analyses 

If an isolated track is allowed to form a jet, both final state, hhqq and hh£ + £~ are characterized 
by six jets, and the hadronic system is characterized by a high b content. 

The analysis works in three steps. The preselection mainly based on event shape variables 
consists in removing most of the two- and four-fermions final-state contributions of the back- 
ground. In a second step the b-content as well as the di-jets reconstructed masses are used 
to reject the three-bosons final-state contribution. Finally, all the relevant informations are 
combined with a multivariable method. 

The preselection consists in rejecting events with a charged multiplicity less than 30. The 
thrust value has to be less than 0.9 and the direction of the thrust with respect to beam axis 
should verify cos(0 t hrust) < 0.9 . In order to remove e + e~— >Z(j), the longitudinal component of 
the total momentum of the event (P* ot ) have to satisfy P* ot < 50 GeV/c. The event is divided 
in two hemispheres with respect to the thrust axis and the invariant mass of each hemisphere 
(M hemi ) must exceed 90 GeV/c 2 . Figure |] illustrates the discrimination between the signal and 
the background processes achieved by this variable. 

The event is clustered with the DURHAM scheme. The variation of the y cut is stopped when 
the transition between six and seven jets occurs and the event is thus forced into a six jets 
topology. To eliminate e + e~— ^(7) with a photon emitted in the sensitive region of the ECAL, 
a fraction of electromagnetic energy with respect to the jet energy less than 80% is requested 
for all the jets. 

The numbers of events expected from background processes at the preselection level are 
reported in Table 0. The contribution from e + e~— >Zee or Wez/ processes is canceled. Around 
48. 10 3 background events are expected with a luminosity of 500 fb~ while ~42 and ~7 
signal events are expected for hhqq and hh£ + £~ final state respectively when mh=120GeV/c 2 
is assumed. Major sources of background are processes with at least one top quark in the 
final state. In fact, 99% of the contribution from the e + e~— >Z(7) process is coming from Z— >tt 
and 80% of the contribution from e + e~— >W + W~ process is e + e~— >Wtb. The remaining part of 
those background processes mainly consists in six fermions final state with two or one b quarks. 
However, even if the background contribution is largely reduced, the remaining part is three 
order of magnitude greater than the signal. 

Thanks to the jet energy flow resolution, a direct use of the reconstructed di-jets masses 
is applied. Among all the di-jets masses in the event, the so-called mass minimizes the 
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Figure 3: Distribution of the invariant mass of each hemisphere (M^ emi ) for (a) signal and (b) back- 
ground processes. 



difference with a mz hypothesis. The b content of the system recoiling to the Z ( J g rec ° l1 ) is an 
excellent variable to reject contribution from background processes. In a first stage, a loose 
cut is applied : at least one b jet among the recoiling jets system (£> rccoil >l) is requested. 
Among the four jets recoiling to the Z (i.e. m^), the di-jets masses m 12 and m 34 are defined 
by minimizing ||m 12 -m 34 || while having both m\i and m 34 larger than 100 GeV/c 2 . The three 



di-jets masses are combined in a simple way in order to form the distance DIST=w((mi2 — 



mh) 2 + {fn 3 4 — mj,) 2 + (m 56 — mz) 2 ), where and mz are the Higgs boson and Z masses. 
The distribution of this variable DIST is displayed on Fig. ^ when a tightened cut on £> rcco11 is 
applied (£ rccoil >2). 

The contribution from the various background sources are reported in Table for the loose 
(£> rccoil >l) and tightened (£> recoil >2) cuts. The signal efficiencies are 43% and 35% respectively. 
Most of the background is coming from e + e~— >Wtb and e + e~— >tt processes while the six 
fermions final state contribution is due to e + e~— >ZZZ with relatively small level with respect 
to the previous contributions. The figure of merit defined as s/vb increases from 0.22 at the 
preselection level up to 3, with an integrated luminosity of 500 fb~ . 

Clearly the combination of the the di-jets masses in the DIST variable is a very crude one. 
Thus, in order to take into account the correlation between the di-jets masses and the b content 
of the event, the variables m.12, m.34, m.56 and £F eco11 are combined with a multivariable method. 
Here a neural network (NNet) approach has been adopted. The NNet output varies from 
-1 (background-like) to +1 (signal-like). Restricted to the positive values, Fig. || displays the 
NNet output from signal and background processes while the contributions from the background 
sources are indicated in Table |2|. The signal efficiency is about 36% and s/vb is ~4.2 when an 
integrated luminosity of 500 fb _1 is assumed. The clear discrimination between the signal and 
the background processes could be exploited to reinforce the figure of merit : as an example, 
with the same integrated luminosity, a cut on the NNet output at 0.5 leads to 31 signal events 
while 27 background events are expected corresponding to a figure of merit of ~ 6. Better 
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Figure 4: Distribution of the variable DIST (defined in text) for (a) signal and (b) background 
processes. 



values could be obtained with thightened cuts and this behavior is fully used in the cross- 
section measurement where a fit of the NNet output distribution is performed as indicated in 
the next Section. 



5 Results 

The cross-section of the e + e~— >hhZ have to be measured, in order to extract A^j^. Such 
a measurement is performed by a likelihood fit of the characteristic variable (DIST, B mc ° l1 , 
NN output ) already used in the analysis. The different selections (indexed as B, C and D in 
Table ^ have been tested. The relative error obtained on the cross-section measurement (Aa/a) 
according to the different sets of 'selection/discriminant variable' are reported in Table |3| for 
integrated luminosity from 500 to 2000 fb _1 . A relative error of 10.3% on the cross-section is 
obtained when NNet selection is used and an integrated luminosity of 2 ab _1 . 

With the same analysis used for mh=120 GeV/c 2 , the relative error on the e + e~^hhZ 
cross-section has been evaluated for higher Higgs boson masses up to 140 GeV/c 2 . The relative 
errors are reported in Table |] as well as the cross-section values and expected number of events 
with an integrated luminosty of 500 fb -1 . Only the 'B/i3 reco11 ' set of selection/variable has been 
considered. Due to the decrease of the cross-section the error increases rapidly up to 17% even 
with an integrated luminosity of 2 ab _1 . Since this result is obtained with a standard selection, 
a dedicated analysis based on retrained NNet which take into account the Higgs boson mass 
would improve the performances. 

To derive the relative error on Aj,Vh (AA/A) the relation between Aj^h and the cross-section 
has to be taken into account according the behavior displayed on Fig. || [|Kj where k is A^j^ 
expressed in unit of Ao, the trilinear Higgs coupling in the standard model. Thus, on this 
Fig. H k—1 stands for the standard model. Furthermore, the incoming beams on that figure 
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Table 2: Numbers of events with £= 500 fb~ expected both for signal and background processes at 
preselection level, standard selections (two set of cut on j3 rec ° l1 ) and multivariable analysis; s/b and 
s/Vb are also indicated. 
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Table 3: Relative error (Act jo) on Ohhz for different selection/variable set and integrated luminosities 

are assumed to be fully polarized so that the displayed cross-section must be reduced by a 
factor two |IJ to be coherent with the present unpolarized study. 

As indicated by the Fig. [I], only one diagram for the e + e~— >hhZ cross-section is sensitive 
to Avu},- The others constitute an irreducible background for the A^h measurement which 
shows up on Fig. |6] as the no n- vanishing cross-section at k — 0. Such an irreducible background 
makes the measurement intrinsically difficult. 

From Fig. [|the relative variations of the cross-section and A^^ around the standard model 
are related by AA/A~1.75x Aa/cr. The relative error on Aj^^ is accordingly reported in Table [5] 
for two sets of select ion/ variable. With an integrated luminosity of 2 ab _1 and NNet output, 
the relative precision AA/A is 18%, and relaxes to 22% with half of that luminosity. 
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Figure 6: e + e -*hhZ cross-section (in fb) as a function of k when the incoming beams are fully 
polarized; m h =120 GeV/c 2 and y/s=500 GeV. 
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Table 4: Relative error (Aa/a) on a hhZ with the 'B/B rccoi1 ' set for different Higgs boson masses and 
integrated luminosities; signal efficiencies ( ehhz ) and cross-sections ( <7hhz ) are reported and (N^hz) 
the expected number of hhZ events with £=500fb~ 1 
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Table 5: Relative error (AA/A) on Xhhh for different selections and integrated luminosities. 

6 Conclusions 

To establish the Higgs mechanism in an unambiguous way, the self-energy potential of the 
Higgs field must be reconstructed. This implies the determination of the trilinear self-coupling. 
The experimental feasibility of the Ahhh measurement has been explored trough a detailed 
analysis of the reconstruction of the double Higgs-strahlung events. It takes advantage of the 
characteristic signature with four b jets and a Z boson, reconstructed either in its leptonic or 
hadronic decay modes. The large four and six fermion background and the tiny signal cross 
section make of this analysis a genuine experimental challenge. An excellent tagging as well as 
reconstruction capabilities of the Tesla detector are then essential. Thanks to a multivariable 
method, based here on a neural network (NNet), a s/Vh=Q could be achieved with 68 signal 
events selected in 1000 fb _1 with a center-of-mass energy of 500 GeV and mh=120 GeV/c 2 ; 
18% of the selected signal events are from hh£ + £~. With help of high integrated luminosity 
(2 ab _1 ), the determination of the e + e _ — >hhZ cross-section with a 10% relative error could be 
achieved leading to a relative error on A h hh of 18%. 
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